Superconducting instability in 3 bands metallic nanotubes 
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Motivated by recent experiments on small radius nanotubes, we study the superconducting in- 
stabilities of cylindrical (5,0) nanotubes. According to band structure calculations, these nanotubes 
possess three bands at the Fermi energy. Using a fermionic renormalization group approach and a 
careful bosonization treatment, we consider the effect of different attractive interactions, mediated 
by phonons, within the Luttinger Liquid framework. We particularly focus on a superconducting 
instability specific to the three bands model we consider for the description of these (5,0) cylindrical 
nanotubes. 
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I. INTRODUCTION 



Superconducting behavior of Carbon nanotubes have 
been observed both in ropes of single walled nanotubesi, 
and in single walledi^ and multiwalledi^ small radius nan- 
otubes grown in a zeolite matrix. The high value of the 
critical temperature in these last two cases raised the 
question of its possible relation with the small radius of 
the nanotubes. Indeed, for this 4 angstrom nanotubes, 
the large curvature induces an hybridization of the a and 
TT orbitals of the carbon atoms^, leading to electron and 
phonon properties different from the larger nanotubes. 
The relation between the high Tc and these peculiari- 
ties have motivated several works, in particular on the 
metallic (5,0) nanotube which constitute the best can- 
didate for the origin of the superconductivitji^. These 
previous approaches include both numerical calculations 
of the band structure and phonon dispersion relationSii, 
and renormalization group approaches either restricted 
to a subspace of the couplings^, or using specific initial 
conditions in the full space of couplings.^ 

In this paper, we identify the different instabili- 
ties of the (5,0) metallic nanotubes in the presence 
of effective electronic attractive couplings mediated by 
phonons. We follow previous approaches on larger 
nanotubea^i"ii^i^i^'^ii4i^'^i^"ii7ii«d3 in using the Lut- 
tinger Liquid framework to describe the low energy be- 
havior of nanotubes. Our approach is based on the band 
structure for cylindrical (5,0) nanotubes provided by var- 
ious methods such as the Local Density Approximation 
(LDA), the GW method, and tight binding calculations, 
consisting in three bands at the Fermi energy^iSfliSLSSiS^ 
Then, we study the perturbations of this band structure 
induced by the residual interactions between the low en- 
ergy fermions. The nature of these interactions is con- 
strained by the specific symmetries of the initial band 
structure, different from usual 3-lcg fermionic ladders 
previously studied in Refs. 2_ 425 .26 27-28 29 3 31 .32 la. 
Using a fermionic renormalization group we identify the 
dominant instability corresponding to each effective at- 
tractive potential. The instabilities we find involve either 
electronic degrees of freedom on a single band, or on a 
two band subsystem, or on the whole three band system. 



We focus on this last case, which corresponds to new 
instabilities specific to the symmetries of the (5,0) nan- 
otubes. This allows to determine the momentum of the 
phonons responsible for the main instability, analogously 
to the proposal for a superconductivity induced by radial 
breather modes in regular 2-bands metallic nanotubes^. 
For the specific instabilities, seen as strong coupling di- 
rections of the renormalization group flow, we use the 
abelian bosonization formalism to identify its nature and 
specify the corresponding dominant correlation function. 
This bosonization description requires a careful treat- 
ment of the so-called Klein factors, a crucial technical 
point in this 3 bands model. We pedagogically present 
this problem and its solution in section Hvl Finally, the 
remaining gapless spin modes are identified using a non- 
abelian bosonization approachi2Si2Si2ii22i22ii£ 

The paper is organized as follows : in section ^ we 
define the fermionic model we consider and the nota- 
tions used throughout the paper. The renormalization 
approach is sketched in the next section IIIII the de- 
tails being postponed to appendix El for readability of 
the manuscript. The strong coupling phases are ana- 
lyzed in sec. IIVI as well as the conventions used for the 
bosonization formalism. The complete bosonized expres- 
sions of all operators and necessary correlations functions 
are given in appendix^ Finally we discuss the validity of 
our results and our main conclusions in section The 
appendix [CI is devoted to some peculiar technical diffi- 
culties of our model associated with the Klein factor of 
abelian bosonization. 



II. MODEL 

A. Band structure of (5,0) nanotube 

The band structure predicted by LDA-DFT calcula- 
tions for (5,0) metallic nanotubes is depicted schemati- 
cally around the Fermi energy Ep in Fig^ It consists 
of three bands^i^fliSiiSSiS^. For a cylindrical nanotube, 
rotational invariance results in the conservation of the 
angular momentum m, and translational invariance in 
the conservation of momentum along the tube k^- The 
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quantum numbers of the three bands are thus determined 
accordingly. In our specific case the two bands with an- 
gular momentum m — ±1 are degenerate and possess 
the same Fermi momentum kp-^, smaller that the Fermi 
momentum kp^ of the band m = 0. Linearizing these 
bands near Ep, we decompose the fcrmion annihilation 
operator into : 

iJoA^, <i>) - e^''-o-^^_o^^(2;) + e-*'=-o-^^_o_^(x) (1) 

and 

+ e+''=-i^-™VL,m.<.(x) (2) 

where 'ipB./L,rn.<T represent the annihilation operator for a 
right (resp. left) moving fermion of angular momentum 
m and spin a. From now on x denotes the coordinate in 
the direction of the nanotube, and </> the angle along the 
circumference (see Fig. . Note that as the Fermi veloc- 
ity vp-^ in the m = ±1 bands in negative, the right moving 
fermions of these bands have a longitudinal momentum 
— /cfi as opposed to the usual situation of the band m — 
where they have momentum +kpg{ see Fig^. We then 
describe the low energy properties of this model by the 
simple Hamiltonian 



E E ^^^-^ 

ff=ta m=0,±l 



(3) 



Note that as the nature of the superconducting instabili- 
ties we discuss in this paper will not depend on the differ- 
ences of Fermi velocities between the bands, we will as- 
sume from now that they are all equals : vp^ = vp„ = vp. 




m=+/-l 



m=0 




L,0 



R,0 



FIG. 1: Schematic representation of the band structure near 
Ep for the 3-bands nanotubes considered in this paper. 



B. Residual interactions 

Then we consider the perturbations around this band 
structure, which can originate either from electronic in- 
teractions not taken into account by the band structure 



calculations, especially in a quasi-lD geometry, or from 
the effective attractive interaction originating from the 
coupling to phonons. Within our effective low-energy ap- 
proach, we consider a minimal model possessing all the 
conservation laws. This results in a perturbative interac- 
tion action which can be written formally as 



a. a' 



where a, b,c,d = 0, ±1 stands for the band indices (an- 
gular momentum). In this expression, as is usual in ID 
systems, the first part correspond to the back scattering 
operators, and the second to forward scattering opera- 
tors. The forward scattering part can be decomposed 
into 54 processes and 52 processes^ii^ The 54 processes 
only renormalize the velocities of the particles^i*^ Since 
we have neglected the velocity differences, consistency 
requires to also neglect the 54 processes. We use the con- 
vention for the action that the partition function of the 
system is written a.s Z = J dipdip''' exp{—S[■^p,1p^), and 
that repulsive (resp. attractive) interactions between the 
fermions correspond to 5*-^\.g^^-' > (resp. < 0). 

To proceed, we must use the symmetry of the problem 
at sake to select out of all the couplings in Q only those 
fulfilling the required conservation laws. As the results of 
band structure calculation suggest, the Fermi wavevector 
/cfi of the two bands ±1 is different (and incommensu- 
rate) from the Fermi wavevector kp^ 7^ kp^i of the band 
with angular momentum m = 0. Interactions must pre- 
serve both rotational invariance and translational invari- 
ance, i.e conserve the total angular momentum m and 
the total momentum k^- To classify these interactions, 
we follow the notations of Refs. |43i and ^ Note that 
whereas this model has some superficial similarity with 
the 3-leg ladder modelSiiSSiSLSSiS, it differs from it by the 
symmetries as all three Fermi momentum are different, 
as opposed to the present case. 



1. Interactions in the band m — subsystem. 

The first two allowed interactions are the usual back- 
scattering and forward scattering interactions in the sin- 
gle band m = 0. The associated fields are denoted by 
g^^^ and g'^^^ and they correspond to the action 




(2) 
1 



g. 



(1) 



(2) 



(1) 







m 





g, 



(2) 



FIG. 2: Formal representation in the kx,m plane of the con- 
sidered interactions in the two band m — ±1 subsystem. 



2. Interactions in the two bands m — ±1 subsystem. 

The next group of interactions we consider are the 
forward and back-scattering couplings in the subsystem 
consisting in the two degenerate bands m — ±1. This 
corresponds exactly to a two leg ladder with degenerate 

bands.'6i44i45,46,47,48,49,5Q,51,52,53,54,55,56 ^ote that even 

this 2 band subsystem differs from the usual description 
of larger nanotubes, which possess right (or left) moving 
fermions at both +kp and —kp.^'^ The interactions are 
depicted schematically in Figure |21 The explicit part of 
the interacting action is 



our conventions for the signs of the coupling, they read 



Ao/±i) 



(7) 



a.a' 

- /^'^ (V'ko,.4,+i,.'V'L,+i,.V'fl,o,.' + (+1 -1) + h.c 
+ " (V'fl,o,aV'Lo,a''/'L,+i,a'V'fl,-i,<T + (+1 ^ -1) + h.c^ 

- ^ {^R, + l,a^i,~l.a'^Lfl.ai^Rfl,a' + ( + 1 ^ "l) + ^-C-) 

All these couplings are depicted schematically in Fig. |3| 
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FIG. 3: Formal representation in the kx,m plane of the con- 
sidered interactions between the m = band and the two 
band m — ±1 subsystem. 



:(±1) 



a, a' 

(fR.. + l.a'l'\,-l.a'i'L.-l,a'^R, + \,a' + 
+ 5f^ (V'fl. + l.^V'I -l,a'V'L, + l,^''/'fl,-l.<T + 

- ('/'k + l,'TV'I,-l,a'V'L,-Hl,^V'fl,_l,<x' + 
+ 52^^ (V'fl, + l,^V'i,-l,^'V'L,-l,<T''0fl, + l,^ + 

- 54^^ (V'fl, + l,^V'L + i^<,'V'L, + l,<x^J?„ + l,<T' + 
(^k + l,^^L + l,<T'^L. + l.<T'V^i?„ + l,<T + 



(6) 
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3. Interactions between the m = band, and the two bands 
m = ±1 



The last group of interactions, specific to the model 
we consider, corresponds to the interactions between the 
band m = and the two band subsystem m = ±1. With 



III. RENORMALIZATION GROUP STUDY 

A. Derivations of the scaling equations 

Having defined explicitly the action describing our 
model, we will now study its low energy behavior us- 
ing the renormalization group formalism. The standard 
procedure for one-dimensional fermionic model is imple- 
mented by using the operator product expansion formal- 
ism (see e.g Ref. l58l chap. 5). The product expansion of 
the four fermions operators appearing in the perturbative 
expansion of the partition function reads formally 



{'^'^R,a^L.b'^LAR,d'lp\^e^L.j'^L,g'lpRM) = 

SbgSde , ,t ,+ , , , 



(8) 
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where we have used mixed labels a,b,c,d,e, f,g,h for the 
band m and spin a. In this expression, Za stands for 
X — iVaT. Within the approximation vp^ — vp, all the 
prefactor will produce the constant 

1 I dzdz dl 

where a is a real space ultraviolet cutoff. Thus, specifying 
the operator product expansion ||SJ| to the interactions of 
our model, we obtain the renormalization group equa- 
tions to second order in the couplings g'f\fi,u, v. They 
are given explicitly in formula lASp in appendix^ They 
only differ from those of Ref.j2,by an extra term 2(u^+f ^) 
in the equation for dig2 ■ 



B. Renormalization flow integration 

The scaling equations ljA3|l admit asymptotic solutions 
of the form 

g»^''^(^)- (//V- +^((^-^*)'"")- (10) 

However, a direct analytical solution in the full param- 
eter space is not tractable. Hence we have numerically 
integrated these equations for different initial values of 
the couplings. Our strategy was to choose some reason- 
able perturbative initial point in the parameter space, 
and to study the instability from this starting point oc- 
curing upon increasing the strength of a given attrac- 
tive interaction. We have done this for all the possi- 
ble attractive interactions, i.e the different phonons cou- 
pling the different electronic branches of the model : 

/^^^ 5p\ 54^^ "^^^ We have also checked 
that the results were independent of the initial perturba- 
tive point chosen. 

Each of these instabilities corresponds to a strong 
coupling direction where at least some of the couplings 
(7p^ , , M, V diverge at a finite scaling length I* . Thus we 
characterize each strong coupling direction by the subset 
of the most diverging couplings, namely those with the 
largest power /i^j in Eq. (|10|) . Indeed, we have found that 
for the considered directions, while the dominant cou- 
plings always diverge with an exponent = 1, there 
exist other couplings diverging with smaller exponents 
Hij < 1. These couplings with weaker divergences are 
expected to give rise to anomalous scalingSiSS but not 
to modify the strong coupling phases. 

The nature of the instability corresponding to a given 
strong coupling fixed point will be identified in the sec- 
tion llVI by bosonizing the model in the subspace consist- 
ing of the dominant diverging couplings. We will focus 
particularly on the instabilities specific to the three band 
model. 



1. Single band or two-band model instabilities 



We first list the instabilities of the band m = 0, and 
two bands m = ±1 subsystem. These instabilities are not 
specific to the presen t mode l, and have been previously 
studied (see e.g. Refs. l6ll63 and references therein). The 
first instability is obtained for a negative gi , i.e an attrac- 
tive interaction in the m — band. The corresponding 
asymptotic fixed point, at which gi and g2 both diverge 
to — oo, is the well known instability of the Luther-Emery 
modeliiiSigi 

Three different instabilities affect only the bands m = 
±1. Upon decreasing g[^'' to negative values, we find that 
the dominant diverging couplings are gi^\ g''^\ g2^\ g^^ 
and g^\ While gl^'' and g^^'' flow towards — cx3, g^{^\g^'' 

(2) 

and gl flow towards -l-oo. The phase associated with 

(2) 

a negative coupling gl ' < is described by : dom- 
inant divergence of g^^ +oo and g^\g4'\g^^ 
—oo. Finally a negative g^^ < induces the phase 
(72^\ g2^\ 34^^ ^4^^ — > —OO. All these strong couphng 
fixed points correspond to the superconducting phase 
of a fermionic two leg ladder, associated with different 

SymmetrieSi^6,44,45,46,47,48,49,51,53 



2. Three bands instabilities 



We now focus on the phases induced by attractive 
interactions specific to the 3 bands nanotube model, 
namely f^^\u and v. In all three cases, the dominant di- 
vergent couplings are 52 , /i , /2 , "* , "* , g!t^ , We have 
identified two pairs of asymptotic directions, one induced 
by negative /'■^^ or u, and the second by v. These two 
pairs of strong coupling directions differ only by the sign 

of the asymptotic v{l*). In both cases, 52, /i, /2, 92 ^ow 

(2) (2) 
towards —00, g\ towards +00 and g^ ±00. The 

first direction corresponds to u — > cx), and the second to 
V —^ —00. Two numerical flow obtained by slowly de- 
creasing either m or to negative values are shown in 
figures 0] and |S1 

It is instructive to analyze further the renormalization 
fiow by focusing in the subspace of the dominant cou- 
plings (72, /i , /2, ffP'' , g^^ , g^^ , V. Indeed this subspace is 
stable under the RG equations (jA3|) . When restricted to 
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(2) (2) (2) 

FIG. 4: Numerical renormalization group flow showing the instabihty induced by a negative u in the Y ~ gl ~ 92 ,9i plane 
(left), g^\v plane (middle), and the ratio u{l)/v{l) as a function of I (right). These flow corresponds to the following initial 
values g'-^'> = 0.17, ff'^^ = 0.3, p'> = 0.1, = 0.09, gf^ = 0.13, = 0-15, g^^^ = 0.05,5^^^ = 0.2, gf' = 0.08, gf^ = 0.1, and 
V = 0.1. u was taken to negative values in step of 0.03 : u = — 0.03i for i = 1, 5. 




FIG. 5: Numerical renormalization group flow showing the instability induced by a negative v in the Y = g]^ —§2 ,§[ plane 
(left), g[^\v plane (middle), and the ratio u{l)/v(l) as a function of I (right). These flow corresponds to the same initial values 
as in Fig. 2] except that u was held to fi = 0.1 and v = — 0.03i for i = 1, 5. 



this subspace, these equations read 







(11a) 






(lib) 






(He) 






(lid) 


a,/«=-2(/«)2 




(He) 


aJ(2) = -(/(^))' 




(llf) 


dii = - (4/(1) - 


2fi^ + ~gf^+gi^- 


^5^)^^ (Ug) 



These equations possess two scaling invariants : C — 
2~gf + 2gf ) - 5(2) and D = 2/(2) _ /(i) + g(2). Let us 
start by considering the RG flow in the subspace v = 
0. Introducing the variable Y = g\ — g2 , the RG 
equations reduce to those of Kosterhtz and Thouless : 

= 2(5(2^)2 ; a,5i'^ = 2g(2V (12a) 
=-2(/(i))2 (12b) 

and g(2) is a flow constant. The asymptotic solutions 
are thus the two directions (A and B on Fig. O Y(X) ~ 
l/(2(r - 0) and ^ = ±l/(2(/* - 0) and the line (C) 
^p-* = 0, coupled to the solutions /(^^ = or /(^^ = 
— 1/(2(/*— /)) . The solutions corresponding to — or 



/(^^ = are easily found to be unstable when introducing 
a small v. 




B 



FIG. 6: Schematic renormalization group flow in the g\ ,Y = 
ffF' ~ 92^ plane, for = 0. 

The scaling behavior of v can be deduced from inspec- 
tion of equations (|ll|l . Except for extremely large initial 
value of 54 ,v will always end up diverging to ±c», in the 
direction given by its initial sign. The scaling equation 
(|llb|l shows that an increasing v(]) leads to an instability 
of the asymptotic direction A. The only remaining scal- 
ing direction driven by v corresponds to the point B in 
Fig. El The other possibility corresponds to strong initial 
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back scattering 


forward scatt. 




J(2) 


(2) 


(2) 
9-2 


(2) 
94 




V 


— 00 


— CO 


— oo 


— CXD 


— OO 


+ 0O 


±oo 



TABLE I: Table of the dominant couplings and their asymp- 
totic directions corresponding to the 3 bands instability on 
which we focus. 



(2) 

intra-band repulsion ^4 , leading to a large and positive 

(2) 

g\ . These strong interaction naturally forbid supercon- 
ducting interactions, and we will not consider them in 
the following. 

To complement the above analysis perturbative in v, 
we have analytically solved the equations Hll|l with the 
ansatz g{l) ~ Ag/{1 — I*) for all couplings. We have 
found that besides the v = Q fixed points, there ex- 
ist three strong couplings directions corresponding to 

~9^'\f^'\f^'\9P - -00, ^ +00 and ~g['^ - ±00 
(together with v — s- ±cx)). The first direction corresponds 

to the direction = l/(2(/* - /)) (point A in Fig. 

(2) 

which is the point induced by a very strong initial g\ 
discussed above. Both the last two directions correspond 
to the same limiting sign of the coupling constants. They 
differ only by the numerical value of the ansatz param- 
eters, and correspond to the main instability discussed 
in the following. We have thus identified a new instabil- 
ity, specific to the 3 bands model we consider. We now 
turn to the bosonization formalism to identify the nature 
of the phase corresponding to this renormalization flow 
direction. 



IV. BOSONIZATION AND NATURE OF THE 
INSTABILITIES 



The purpose of this section is to identify the nature 
of the instability for the previously identified strong cou- 
pling direction listed in table This will be achieved us- 
ing the bosonization formalism within the subspace cor- 
responding to the dominant couplings. We will pay spe- 
cial attention to the proper definition of so-called Klein 
factor. We first start by defining our conventions on the 
non-interacting three band model defined in section UTI 



A. The "Condensed Matter" Bosonization 
Dictionary 

In the standard "condensed matter" bosonization pro- 
cedure, we represent the annihilation operators of right 
and left moving fermions, defined in aaSS^SSiSi: 




-i*H,m,<T(a:) 



(13a) 
(13b) 



where we introduced Majorana fermion operators (the 
so-called Klein factors) ?7fl/L,m,CT that satisfy: 



m,o" 1 ^i?,,m' ,0"' }+ '^^rn.m' ^a,cr' (14r) 
m,cr ; ^L,m',o"' }+ ^ ^^^771,771' ^ a. <t' (14b) 



'}+-o. 



(14c) 



Note that in this convention, we introduce one Klein fac- 
tor per set of quantum numbers (±kp, m, a). With these 
anticommutation relations, the proper anticommutation 
relations for the fermion operators defined in (|13|l are 
satisfied with the following commutation relation of the 



fields <i> 



R/L 



[^R,m,a{x), ^R,m'.cr'{x')] = inSm^m' S^y, a' sign{x - x') 
[^L,Tn,a{x),^L,rn'.cr'{x')] = -i'!TS,n,Tn'Sa,a'Sign{x - x') 

')]=0. (15) 

With these conventions, the bosonized non-interacting 
Hamiltonian readsMigiaSLSl: 



H 



E 

n = 0,±l 



/ -^^F [{^^R..^,af + (V$L,™.,)'] , (16) 



Finally, the densities of right moving and left moving 
fermions read respectivelySiSSiSl: 



PR,m,a 
PR,m,<7 



R.m.a 



2tt 



2tt 



(17a) 
(17b) 



Keeping only the 52 processes defined in section III Bl 
we express the forward scattering part of the interactions 
in terms of the above densities ()17fl : 

-ff forward =ffF' '^2iPR,l,(TPL,l,a' + PR,-l,a PL -l,a') 
(T.cr' 

+/^^^ 'YyyPRS,'y + PR,-i,>r)pL,o,<T' 

a, a' 

+ {PL,l,a + PL.-l,a)PR,Q.a'] 
+52^^ '^{PR,l,<TPL-l,a' + PR,~l,aPL,l,a') 



+9'''^^^PR,0..>7PL,0,a' 



(18) 



Using (|17|l . these expressions can be reduced to quadratic 

expressions in the fields ^R/L,m,cr- Note that we will 
(2) 

treat the gl ' term below with the backscattering part 
of the Hamiltonian : although it appears as a forward 
scattering term, it cannot be reduced to a density-density 
coupling and its treatment closely follows the one for the 
backscattering couplings. 
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B. Backscattering Interactions 

1. The Klein factors problem 

Whereas the bosonized forward scattering part of the 
Haniiltonian is function solely of the densities, and thus 
does not depend on the convention chosen for the Klein 
factor H14|) and fields (|15|l , the situation is different for the 
back-scattering part of the action 10}. Quite generally, 
these back-scattering operators can be written as 

where we assumed the transverse momentum conserva- 
tion mi -I- TO2 = ma -I- m4 and we defined the product of 
Majorana fermions 

and the product of vertex operators 

(21) 

The usual strategy is to find a representation such that 
the operators 0^m},<j,a' commute with each other, and 
similarly for the products of four Majorana fermions 
J^{m},a,(y' ■ It IS imp ortant to note that the only case 
discussed in Ref. is the one in which all the vertex 
operators 0{m}.(j,(T' are already commuting so that no 
redefinition of the fields is necessary. When determining 
the ground state, it is then possible to consider the Sine- 
Gordon form of the Hamiltonian, obtained by replacing 
the operators M{m},(j,a' by their eigenvalues. 

Obviously, while the 4 fermions operators (|19|l always 
commute with each other, the above condition of inde- 
pendent commutation of the N{m},<T,a' and 0[m},a,i7' be- 
comes more and more difficult to fulfill with an increasing 
number of fermion species. This is particularly true for 
our 3 bands model, corresponding to 12 fermionic species 
(2 spins and 6 Fermi points), and we can check on the 
bosonized expressions derived in App.m Eqs. ||B2|) 
and IjBSp . that this condition cannot be satisfied for the 
operators we consider within the convention defined by 
Eqs. p 3114115(1 . Indeed, the products of four Majorana 
fermion operators are commuting when they have an even 
number of Majorana fermions in common and anticom- 
muting otherwise. In the second case, which occurs for 
our model (see e.g the operator v), the corresponding 
operators 0{m}.CT.(T' also contains an odd number of ver- 
tex operators in common. Since the vertex operators 
associated with fermions are anticommuting when they 
correspond to the same fermion species, we recover in 
the bosonization formalism the commutation of the four 
Fermi operators, but the independent commutations of 
the C{m},CT,cr' and A/{m},cr,CT' is not possible. We thus need 
to change our bosonization convention for this particular 
model. 



2. "Field Theory" Convention 

The problem we have to deal with is thus whether it 
is possible to redefine the Majorana fermion operators 
(|14|l and the commutation relations of the fields H15() 
that appear in the vertex operators in such a way that 
all the new products of four Majorana fermion opera- 
tors are commuting with each other and simultaneously 
all the new products of vertex operators are also com- 
muting with each other. Another possible convention, 
different from the above "condensed matter" convention 
(left /right moving chiral field of the same band commut- 
ing with each other, and one Majorana fermion per Fermi 
point^^), consists in what we will call the "quantum field 
theory" convention. We now introduce a single Majorana 
fermion for a pair of right and left fermion which need 
not have the same quantum numbers (same specie). Cor- 
respondingly, the chiral fields for this pair have a nonzero 
commutator. In this representation, the fermion opera- 
tors are now expressed asS^ 

V'i?.,™,.(a;) = -^e-^*«— (22a) 
i^L.Pin.)A^) = ^^e^*"--^"^?7m,^, (22b) 

where P is a permutation of the band indices (fermion 
species). Now the field commutations relations are mod- 
ified into : 

[^R.m.aix), ^R.rn',cr'{x')] = iTr5m,m'Sa.a'Sign{x - x'), 
[$L,m,<T(x),$L,m',o-'(a;')] = -iTr6,n,m' Sa^a' sign{x - x') , 
[^R,m,cr{x),<^L,m'.cr'ix')] = iTTSm,m'Sa.a' (23) 

and the Majorana fermion operators 77m, o- satisfy: 

{Vm^a : Vni' ,a'} — '^^7n.7n'^{T,a' • (24) 

We discuss the equivalence of these two representations 
in the appendix [CI Let us now apply this convention 
to the present model. We have found that the suitable 
(necessary) permutation P of band indices in Eq. H22|l 
simply permutes the bands -1-1 and —1. 

The remaining interactions f^^\v are conve- 
niently expressed in terms of the following non chiral 
fieldsi^: 

4>m,a = -{^L,m,a + ^R,m,a)- (25) 

These fields satisfy [4>m,ij{x).,4'm' ,(t'{x')] — 

[Om.Ax).d,n'M'{x')] = and [(j)m,a{x),Om' ,c,'{x')\ = 

*f '5m,m"5CT.(T'Sign(x' — x). Taking the derivative with 
respect to x' and introducing Iirn.a{x) = ^dx6„i^a one 
finds [4>m,a{x),Ii.m,a{x')] = i5{x — x'), showing that the 
fields Tim, a and 4>m,a are canonically conjugate. It is 
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convenient to introduce the total and relative "charge" 
and "spin" fieldsi^: 



7m, T + Om.i ) 



(26a) 
(26b) 
(26c) 
(26d) 



Finally, the following rotation for the charge modes will 
simplify the expressions of the interactions : 









1 






V3 






1 


1 










v 


i>c,B j 


\ 76 


1 



1 

\/3 





(27) 



The same rotation is performed for Od so that canonical 



commutation relations are preserved. A different rota- 
tion must be performed for the spin modes: 



1 


1 


\/2 




1 


1 


V2 


V2 



(28) 



With these notations and within the above field-theoretic 

(2) 

representation, the bosonized expression of the gl term 
becomes: 



(2) 

9i 



(tto) 



2 , dx cos 26c, A ( cos 2ds 



+ cos2(f>s-r)i.Am,lV-i.lV-iA ]■ (29) 



After some simple algebra, we find for the /^^^ part of 
the Hamiltonian the simplified expression 



/(I) 



(27ra)2 

and for the only remaining v coupling, the bosonization expressions reads 
2v 



, (30) 



(tto) 



2 S VilVilVoiVo^ 



cos cos0cA cos(/)s_ cos(\/20sO + (t)s+) + sm^ViOcB) sinOcA sin0s_ sin(\/20so + 0s+) 

(31) 



cos(\/30cs) cos^cA cos^s- cos(\/26'so ~ ()s+) — sin(\/30cs) sin^cyi sin6's- sin(\/26'so ~ 

I 



We immediately observe that the change of bosoniza- 
tion convention results in an important simplification. 
First, in this field-theoretic representation, the Majo- 
rana fcrmion product in the gi term is now commuting 
with the Majorana fermion products that appear in the 
/^^^ and the v terms. This allows for a simultaneous di- 
agonalization of all these remaining Majorana fermions 
products. Representing the 2 Majorana fermions prod- 
ucts as pseudo-spins : jyi.i^^i.x — ifif , ?7-i.|?7-i.x = id'^i, 



Vo,]Vo,i = ^^■g, and choosing the -1-1 eigenvalues of the 
CT^, we obtain the final bosonized action, which takes 

the Sine-Gordon like form. The first g^ term reads 



-p—^ [ dx cos 29^ A (cos 29 + cos 2(j)s- ] , (32) 
the ' term simplifies into 



p(i) 



(7ra)2 



<?!>s,+) - V2{9sfi + 0s,o)l cos(6is,- " 0s,-) \, (33) 
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and finally the v term can be written as 

2v ( 



cos 0cA COS (f>s- cos(-\/2(/)so + (t's+ ) + sm{Vi0cB ) sin 9cA sin 0s- sin(A/20sO + 4>s+ ) 



cos(V30cs) cos(?cA cos cos{V29so — 9s+) — sin(\/30cs) siiiOcA sin 0s- sin{V20so — 0s+) \ (34) 



r 



Analysis of the strong coupling fixed points : 
superconducting instability 



Having obtained the above expressions Eqs. 



we are now ready to characterize the phases correspond- 
ing to the strong coupling directions identified in the 
renormalization study of Sec. IIII B 21 We will only fo- 
cus on the new instability specific to a three band model. 

(2) 

This instability corresponds to a divergence of 3^ —> 
—00, /^^^ — > -l-oo and v — > ±00 (see table IJ). The two 
signs of V are possible, depending on the driving attrac- 
tive perturbations {e.g u or v). From (|32|l . we find that 

(2) 

large negative values of g-^ induce a locking of the field 
0c, A — 0. 0s- and (f)s- being dual to each other, no fur- 
ther information on the spin part can be gained at this 
point. The /^^^ being a pure spin-current interaction 
part, we will postpone its analysis to the next section. 
And finally, plugging the result 0c, a — into 134|) . we 
find that large values of v will induce a locking of the 
charge field 6c. b to 



^ if w > 

v3 

if w < 



(35) 



Thus we find that in the ground state corresponding to 
this instability, the fields (pc,A and 0c, b are locked so as to 
minimize the condensation energy. This implies that the 
corresponding charge degrees of freedom develop a gap. 
The total charge remains gapless as a result of the global 
U(l) symmetry. Thus, only a single charge degree of free- 
dom remains. The analysis of the spin degrees of freedom 
is more difficult, and is done in the next section. How- 
ever, knowing which charge modes are gapped already 
enables us to determine some of the order parameters, 
and find the corresponding nature of the instability. 

In any case, the long range ordering of the charge 
fields 0c. B and 4>c,A has important consequences. In- 
deed, it is seen from the bosonized expressions (IIjl5lij21|) 
of the charge density wave operators that none of them 
can develop quasi-long range order. On the other hand, 
superconducting fluctuations are strongly reinforced by 
the ordering of the charge fields, as can be seen on the 
corresponding expressions ljB7|l derived in the appendix 
IbI Hence we have analyzed the new instability as being 
driven by superconducting fluctuations. 

It is worthwhile to contrast our results with those pre- 
viously obtained on three-leg ladders.-^ In our notations, 



it was found that in the three leg ladder system, the only 
charge field developing long range order was 0c, a- Here, 
by contrast, we find that two charge fields are develop- 
ing a long range order, 0c, a and 9c. b- The difference 
of behavior of the three-band nanotube and the three- 
leg ladder is a consequence of the equality of the Fermi 
wavevectors of the band of angular momentum ±1 which 
is itself a consequence of the rotational symmetry of the 

tube. This equality of wavevector allows extra interac- 

(2) 

tions between the bands ±1 such as gl and between the 
two bands ±1 and the band (such as u or v). The exis- 
tence of these interactions is driving the system to a new 
fixed point. As one more charge modes is gapped in the 
nanotube compared with the three leg ladder, the rein- 
forcement of superconducting fluctuations is expected to 
be a stronger effect in the nanotube. 



D. Effective low energy spin theory for the 
instabilities 



The charge modes of the nanotubes being gapped, 
apart from the global decoupled charge mode, the corre- 
sponding low energy description of the instability we con- 
sider consists only of the spin modes. Further progress 
in the understanding of this theory can be made by in- 
troducing the pseudo-fermion creation and annihilation 
operators: 




V27ra 



(36) 



and the associated Majorana fermion operators (v 
R,L): 



V2 
1 



(-Cl - Kv,2), 

(Ca + Ku,o)- 



(37) 
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(2) 

The interaction term proportional to is then rewrit- 
ten as: 



„(2) r 

2i— / da; cos 26'c,ACfl,oCL.o 
TTa J 



(38) 



(2) 

we see that the interaction term proportional to g\ gives 
a non-zero mass to the C-R/l,o Majorana fermions, while 
leaving the C,r/l,?, fermions massless. This Ising critical- 
ity is a consequence of the self-dual^SiSS character of the 
interaction (|32|l . In the context of two-leg ladders, these 
self-dual interactions have been discussed in Refs. liglTOl 
Moreover, the interaction term proportional to f^^^ 
can also be reexpressed in terms of the fermion fields 
C-R/L,i.2,3- Indeed, we have the relations: 



COS{0 s- + (f)s-) = -^(Cl, 2Cl,3 + <L,3Cl,i): 
■ COs(e's- -<t>s-) ^ -i(Ci?.,2Cfl.3 + «Ci?.,3Ci?.,l)- 



(39) 

In fact, this representation is well known2i and has been 
used to study the two-leg spin ladderi^SiS and the two- 
channel Kondo effectSLZi. Using the equivalence be- 
tween Majorana fermions and the two dimensional Ising 
model, it is also possible to reexpress the interaction v us- 
ing order and disorder parameters of the quantum Ising 
modelS'£S'££'£2*£^. Indeed, one has the relationsSli: 

cos(/)s+ = ^il^i2, 

COS0S+ = inicrifi2, 

sm9s+ = -«K2^ia-2, 

sin(^s+ = -iKiK2CTicr2, (40) 

and similar relations for and 9s- with (/ii,cri) — > 
(MSiCs) and (/U2,cr2) — > (moiCTq)- With these relations, 
we easily find that: 



cos (ps_e ^■ 
cos Os-e''^' 

,(2) 



M3Mo[MiM2 + KlK20-l(T2] 
«K30-3Mo[«K 10-1^2 + ^J■lK2(T2] (41) 



Noting that g\^' — > — oo in Eq. (|38|l implies that /io de- 
velops long range order for (Oc^a) — 0, we find that in the 
low energy limit the expressions in (|41() reduce to: 



cos q. 
cost 



i4's + 



M3[MlAi2 + KiK2(Ti(T2], 

ZKgCTg [iKlCri^2 + Ml'«2CT2] (42) 



Introducing the Pauli spin matrices Tc = — -^^abci^bi^c 
these expressions are easily seen to reduce to the expres- 
sion of the spin 1/2 primary fields of the SU {2)2 Wess- 
Zumino-Novikov-Witten mode l?^i?Ti??i7ii7? in terms of 
Ising fields. Moreover, the expression of the spin cur- 
rents also reduces to the SU{2)2 form. 

Thus the theory describing the spin excitations at low 
energy reduces to a SU{2)i WZNW model (that de- 
scribes the spin excitations of the band 0) coupled with 



a SU{2)2 WZNW model (that describes the spin excita- 
tions of the bands ±1) by a term: 



A / dx ir{gia){x) ■ ir{g2(7){x) 



(43) 



and a marginal current-current interaction term. Power 
counting shows that the term H43|) is relevant with RG 
dimension 5/4. Therefore, it is reasonable to treat first 
this relevant term, as was done in the case of two-leg spin 
laddersiS For analyzing the effect of the interaction H43() 
on the spin spectrum, is convenient to introduce a coset 
representationSi: SU{2)i x SU{2)2 - SU{2)^ x TIM 
where TIM stands for the tricritical Ising model^^. With 
the coset decomposition, we can rewrite the WZNW 
fields as: 



92 = CTTIAigs, 



(44) 



where gk is the spin 1/2 SU{2)k field, ctim is the energy 
operator of the Tricritical Ising Model of dimension jq, 
axiM is the spin operator of the tricritical Ising model 
of dimension ^ . Using the Operator Product Expansion 
of the TIM tTiM<^TiM ~ (<T + <t')tim from Ref. ^2, (p. 
224) we can rewrite the interaction H43(l as: 



A' / (ix(TT/Mtr(.g3cr)(a;) • tr(g3cr)(x). 



(45) 



where only the most relevant term has been kept. The 
interaction is now brought to the form of a self-coupling 
for the SU{2)3 WZNW model. Now, we simplify this 
self coupling by using a second coset^ representation, 
SU{2)3 ~ C/(l) X Z3 where Z3 represents the critical 
3-state Potts model (or equivalently the 3-state clock 
model)'SS. and U{1) represents a free bosonic field de- 
scribed by the Hamiltonian: 



H 



dx 
2^ 



(46) 



where K = 1. This coset representation was used in 
Ref. |8J in a study of the Haldane gap in spin-S chains, 
and from now on our treatment follows this work closely. 
We write the components of the fundamental field as 

with m,m' ^ ±1/2. We have from Refs. I7ll83l 
the relations: 



53 



2 ' 2 J 



^ (Jl, 



(47) 
(48) 



where (Ti,/ii are the order and disorder parameters of 
the 3-state clock model. With this, we can rewrite the 
interaction as: 



A / dxa- 



TIM 



^/2 
PI 



iHfAn (49) 
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Then we use the properties of the 3-state clock modelS: 
(t| = a2 = (72 and similarly with n <-> a, together with 
the Operator Product Expansion aia2 ''^ £■ This allows 
to reduce the above interaction term to: 

a" J dxGTiM [ez, + e-2'\/f-*^2 + e-^'^^^^i^ . (50) 

Now let us make the assumption that the TIM develops a 
long range order and so does the 3-state clock model. Let 
us assume further that the 3-state clock model is in the 
low temperature phase, with disordered. Only the 
bosonic field 4> can a priori remain gapless. In order to 
determine whether (p indeed remains gapless, we have to 
consider the perturbations generated by the disordered 
operators /ii,2- It is straightforward to see that these 
terms yield a perturbation: 

AoCOs2\/60 (51) 

for the U{1) theory. The operators of the SU{2)^ theory 
then reduce to^: 

n+ = tr{ga+) - e'^V^^ (52) 

= tr{ga') - {e-'VI-P^^ + e'VfV2) (53) 

X (e-2»\/iV2 + e-2^^Vi) ~ e-'^^ + e*^^ 

The SU{2) symmetry of the system imposes that n+ and 
have the same scaling dimension. Therefore, at this 
new fixed point, one must have K = 1/3. Hence, after a 
rescaling <j) — ^/\/3 and 9 — V36, K — 3K the expres- 
sions (|46|) and (|52|l reduce to the ones of the SU{2)i 
case^i, with a perturbation (|51|l which is marginal. 
Two regimes are possible depending on whether Aq is 
marginally relevant or marginally irrelevant. In the first 
case, a spin gap is obtained. In the second case, no 
spin gap is obtained and the system has the same spin 
correlation as a free SU{2)i model up to logarithmic 
correctionsi^i*2Si2& In order to predict which phase is re- 
alized, we have to consider the flow of the Luttinger ex- 
ponent of 4>. If the fixed point is approached from the 
side where the perturbation H51|) is irrelevant, then we 
can expect the fixed point to be stable. Since at the ori- 
gin the Luttinger exponent is K = 1, and at the fixed 
point it is iiT = 1/3, the flow is indeed on the side where 
(|51|l is marginally irrelevant. Thus, we find gapless spin 
modes at the fixed point with both triplet and singlet 
superconductivity. However, the logarithmic corrections 
induced by the marginally irrelevant perturbation at the 
SU{2)i fixed point are known4i4Si2& to lead to domi- 
nant triplet superconductivity fluctuations. Triplet su- 
perconductivity is thus naturally expected in the present 
case. In this respect, we note that a similar situation 
arises in single chain system where the renormalization 
group predicts dominant triplet superconducting fluctua- 
tions in the vicinity of a spin density wave phase. Since 
in Refs. 2 3 the superconductivity appears to be sensitive 
to the application of a magnetic field, it is likely that the 



intertube coupling tends to better stabilize the singlet 
superconductivity with respect to the triplet one. If we 
assume that /ii is ordered, then we find that (j) is also 
long range ordered. As a result, n+ is short range or- 
dered while or e = tr[g3) is long range ordered. Since 
the system has to be rotationally symmetric, the only 
solution is to have e long range ordered and short 
range ordered. In that case, the system has a spin gap, 
and only the singlet superconducting order parameters 
exhibits quasi-long range order. 



V. CONCLUSION 

In conclusion, by means of fermionic renormalization, 
abelian and non-abelian bosonization we have analyzed 
the low energy properties of a three band one dimensional 
model deduced from the band structure of cylindrical 
small radius (5,0) nanotubes. We have found that this 
system possess a specific instability, besides the usual sin- 
gle band and two band models instabilities. This insta- 
bility corresponds to the development of superconducting 
fiuctuations in the nanotube. Within our approach, in 
the absence of a spin gap, triplet superconductivity fluc- 
tuations are expected to be dominant du to logarithmic 
corrections, with subdominant singlet superconductivity 
fiuctuations. 

This new instability is tightly related to the symme- 
try of our three band model, and more precisely to the 
new couplings u and v. In our model, in the pres- 
ence of these couplings, the spin excitations are either 
fully gapped leaving only a CISO phase as in the two-leg 
ladder^6,44,45,46,47,48,49,50,5i,53 ^j. ^hey are described by a 

SU{2)i WZNW model leading to a ClSl phase as in a 
single chain Hubbard model. This is in contrast to previ- 
ous studies of three legs model with different symmetries, 
which included only 2 bands couplings as opposed to the 
3 bands couplings u or v : in these models, a C2S1 phase 
was found^^'^'^'^^. Technically, this difference lies in the 
ordering of the field 9c. b, directly related to the presence 
of the V coupling. Note that the 3 bands nature of the v 
coupling also induced the technical problem of the Klein 
factor discussed in this paper. 

Let us finally relate our results to previous studies on 
the (5,0) nanotubes. In Ref. 0, only a subset of the cou- 
plings of the present model was considered, which did 
not include the u and v term. Hence this new supercon- 
ducting instability was not discussed. In Ref.0, Gonzalez 
and Perfetto studied the same model as ours, by means 
of a renormalization group procedure. The nature of the 
phase was determined via the scaling behavior of correla- 
tions functions, as opposed to the bosonization procedure 
used in this paper. In Ref. 9, it was found that the domi- 
nant instability would be a charge density wave coupling 
the bands ±1, with subdominant spin density wave fluc- 
tuations, whereas we find that charge density wave fluc- 
tuations are suppressed. The origin of this discrepancy is 
that in Ref. 9, only speciflc initial conditions were consid- 
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ered, with initial values of some couplings so large that 
they render a one-loop renormalization group approach 
questionable. It appears likely that the initial conditions 
chosen in Ref. strongly favor a two band instability be- 
tween the band ±1. Indeed, a divergent g'^^^ as we found 
in Sec. IIII BTI indeed leads to a reinforcement of charge 
density wave fluctuations between the bands ±1. 

Along these lines, let us mention that it is difficult 
to determine which of the possible mechanisms, includ- 
ing the one proposed in this paper, actually takes place 
in a (5,0) nanotubc. Indeed, as opposed to theoretical 
approaches of conventional larger nanotubes, our one- 
dimensional electronic model is using the band struc- 
ture provided by ab-initio calculations as an input. Since 
these methods already include a renormalization of the 
band structure by a fraction of the electronic interactions, 
an estimate of the bare coupling in our one dimensional 
model based on an unrenormalized Coulomb interaction 
as in Ref. is likely to lead to misleading results by 
overestimating the effect of some interactions. As a re- 
sult, we can only propose a classification of the various 
fixed points at weak coupling and characterize the possi- 
ble scenarii, with the usual hypothesis in one-dimensional 
systems that the weak coupling behavior and the strong 
coupling behavior are continuously connected^^ A re- 
lated remark is that the gaps calculated in any weak 
coupling approximation are generally very small. How- 
ever, in the real system, where interaction strength can 
be expected to be comparable to the bandwidth, since 
Luttinger exponents are found in the rangeSi 0.2 — 0.5, 
the real gaps can be much higher than those estimated in 
a weak coupling treatment. Thus, the present treatment 
cannot lead to a realistic estimate of critical tempera- 
tures. 

Another aspect of the physics to consider is the pos- 
sibility of a pseudo-Peierls transitions in this small ra- 
dius nanotubes^. Indeed, the approach of this paper 
is based on the band structure of numerical approach 
which did not consider the possibility of a cylindrical ge- 
ometry breaking. If such a phenomena was to happen, 
an indication of strong electron-phonon coupling in the 
system, other mechanisms for superconductivity could 
occur, but their description is beyond the scope of the 
present paper. Finally, let us mention that the experi- 
mental results on the superconducting transition in these 
small nanotubes^'^ suggest that a real three dimensional 
superconducting phase transition takes place. A com- 
plete understanding of these results must also include a 
coupling between the nanotubes to stabilize the super- 
conducting fluctuations at nonzero temperature. How- 
ever, since the gap in the zeolite matrix is of order 4eV, 
an intertube Josephson coupling term would be a pri- 
ori strongly suppressed by the presence of the insulator 
between the tubes. 

We thank X. Blase for very stimulating discussions 
which initiated this work. P. Pujol is also acknowledge 
for insightful remarks on the non-abelian bosonization 
approach of section llVDI 



APPENDIX A: DERIVATION OF THE RG 
EQUATIONS 

To express the RG equations in a simpler form, we first 
define rescaled couplings as 



9^ 



1 .U) . /W^_L/W (Al) 



2nvp 



2itvf ' 



1 1 
u = u : V = V 



(A2) 



2TrvF 2'itvf 
The scaling equations read in terms of this couplings 

=-2(g(i))2-4w (A3a) 

di~g^'''^ =-{g^^^f -2u^ -2v'' (A3b) 

a^ffW = - 2[g^l^f 2gf^~g'i^ - 2™ (A3c) 

dra?' ^ - 2~g['^gi'^ - 25f + 2.9f 'gf ' - ^ 

(A3d) 



^ - 25^3^ - 25^5^ -f 2gr 5^ 

- 45«g« + 25«5f ) - 2ui (A3e) 

dfgi'^ - Cgi'^)' - (fff V - (A'^r ~u'-i' (A3f) 

dw'P ^ + 2~gf^~g^^ ~2{g^^yf ^2{g^^f (A3g) 

dfaTHg'^'f-ig'Pf (A3h) 

(A3i) 
(A3j) 



~(i);,(2) o;;(i);;(2) , o;i(2);;(i) 



=-2(/(i')2 + 2u5-2{)2 
d^u = [2J^^-~gf^-g^^-~g^^^) 



(A3k) 



div = - (^-2/^^^ + g^-' -h ffi ' -h 52 ' 1 " 

- (4/(1) - 2/(2) + .#) + 5(^) + .gf ) V (A31) 

APPENDIX B: BOSONIZATION 

In this appendix, we provide some technical details on 
our bosonization approach. We first give the bosonized 
expressions for the three interband operators f^^\g^^ 
and V within the usual "condensed matter" convention. 
These expressions show why this common convention is 
not suitable for the present analysis. Second, we provide 
the detailed expressions of the various order parameters 
of our three band model within the field theoretical con- 
vention used and defined in the text. We use the defini- 
tion and conventions of the section Hvl of the paper. 



1. Derivation of the bosonized form of interband 
interactions 

We consider the interactions f^^\g^'^ that involve a 
transfer of fermions between two different bands. These 
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terms can be read from the action We use the con- densed matter convention of section HV Al 
I 

(2) 

The gl term can be expressed in the form : 



(2) 

9i 



"■^^ [e''^'^'*'"-~VRA,<yVLAMVL.-lMVR.-iM + e''^'^^'--T]R.i,„r]L,i.-aVL,-l-aVR.-i,'y] + H. c. 



(27ra)2 

The contribution of /*^^) is composed of a term quadratic in the $ fields and a term: 



(Bl) 



(27ra)- 



(B2) 



Note that in nonabehan bosonizationSSi^SiiiiS the term /^^^ is seen as the interaction of a SU{2)i current associated 
with the spin excitations of the band with a SU (2)2 current composed of the sum of the SU {2)i currents associated 
with the spin excitations of the bands ±1. And finally, the contribution of v reads: 



(27ra) 



+ e'^^^^''°+'^^-^-''^'-^VRA,aVL,-l,^aVR,0^^aVL^O,a 



_^^ — iV39a,B^ — i4>i:,A 



+ e^-(^'f>^■0+<P.,^+^9.,-)rJLXam.-l,-.VL,0,-.VR,0,a 
I 



H. c. 



(B3) 



2. Order parameters 



perconductivity into bosonized variables : 



In this last section, we switch back to the field theo- 
retical convention for bosonization, defined in the section 
ITVB21 



a. Superconductivity 

We consider the following order parameters the forma- 
tion of singlet superconductivity in the nanotube : 

Oo{x) = "^IpRfi^alpLfi-a (B4) 
a 

Olix) = ^1pRA,ai^L,-l-a (B5) 
a 

0-l{x) = ^ V^_R,-l,cr^L,l,-cr (B6) 



Using the bosonization decomposition introduced in 
the text, we can express these order parameters for su- 



Oo[x) — — e L V J 



^] sinV2(/>s^o 



Oi{x) = — e 



\/3 \/2 V6 



sm(0s,+ 



Oi(x) = — e'^Hw-^+^j sin( 
TTa 



(B7a) 

,-) 

(B7b) 



(B7c) 



Note that we have not considered the triplet supercon- 
ductivity order parameters. Indeed, they naturally pos- 
sess the same charge part than the singlet superconduc- 
tivity operators, and only differ by their spin part. Since 
the spin part is more conveniently treated within the non- 
abelian bosonization formalism (see Sec. lIVD|l . it is not 
necessary to give an explicit expression of the triplet op- 
erators here, since thay can be obtained from the expres- 
sion of the operators (jBTp in nonabelian bosonization by 
the substitution tr{g) — > tr{ga). 
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b. Charge density waves 

Besides superconductivity, one can also expect to ob- 
serve charge density wave order to develop at low temper- 
ature in a quasi- ID system. The various charge density 
wave order operators, labelled by their ordering momen- 
tum vector, are defined as ; 



We have considered only the charge density wave op- 
erators of the form i^li „i ^'4^L,m' ,cr, as those of the form 

"^R m cr'^R,m',a are Current that cannot develop quasi- long 
range order. 





= n ^^L <J 
a 


(B8) 




a 


(B9) 




\ \ " /t / 

a 


(BIO) 


0(kFg-kF^,Ky)iX 




(Bll) 


0(kFg-kF^-Ky){X 


a 


(B12) 


0(2kF^ ,2Ky {X 




(B13) 


0(2kFi,-2Ky{x 


a 




0{2kFo,0){x) = 


— e L J cos 
na 




O{2kF,.0){x) = 


27ra 




O{2kF_,,0){x) = 


27ra 


a 



The corresponding bosonized expressions of these 
charge density wave order parameters is easily obtained 
and read 



0(-kFo-kF^.Ky){x) 



2'Ka 



0,.kF,,.kF,..Ky){x) = ^e'l^^^^^-^^-^+^W.;^-^] Y^-\-hio...^^^ 



i2kF,,2Ky) - 1^ 



0(2kFi-2Ky 



2'Ka 



n/2 



(B15) 
(B16) 

(B17) 
(B18) 



^ J f?o,<T»7i,o- (B19) 
(B20) 
(B21) 



APPENDIX C: EQUIVALENCE OF 
BOSONIZATION CONVENTIONS 



tors used in condensed matter literature amounts to 

writo'6,65,67,89. 



In this appendix, we discuss a general case of ID 
fermions with N "flavors" since the results are of more 
general applicability than the nanotube with three bands 
at the Fermi level that we have considered in this pa- 
per. The bosonized representation of fermion opera- 



1 



'2'Ka 
1 

I 2k a 



(CI) 
(C2) 
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where: 



and one has: 



[*fl,n(a;),*i?,«'(x')] 



i7r(5„,„'Sign(x - a;'), (C3) 
-i7r(5„^„/sign(x - x'), (C4) 
0, (C5) 



and the Majorana fermion operators yy^.n satisfy: 

In the field theoretical literature, an apparently different 
bosonized representation is used^SiSfliSi: 



i'R.,n = 



1 



v27ra 
1 



„-4*H.„(2;) 

^ 'In ; 



i*t,„(x) 



(C7) 
(C8) 



where this time: 

[$it,n(a;),$i?,„'(x')] = i7r(5„,„'Sign(a; - x'), (C9) 
[*L,n(a;),$L.n'(a;')] = -i7r(5„,„/sign(a: - x'), (CIO) 
[^fl,n(a;), $L,«' (a;')] = iTr6n,n' (Cll) 
and the Majorana fermion operators satisfy: 

{lln,Vn'} = 26n,n'. (C12) 

While in the condensed-matter bosonization for a given 
flavor the right and left bosonic fields are made commut- 
ing and there is one Majorana fermion associated with 
the left mover and another Majorana fermion associated 
with the right mover, in the field theoretical representa- 
tion, there is only one Majorana fermion for each flavor. 
The price to pay for this is to make the commutator of 
the chiral fields non-zero. An application of this represen- 
tation in condensed matter physics is the the derivation 
of the bosonized form of the doubled Ising model^SiSLi^. 

To show that these two representations are in fact 
equivalent, let us introduce the conjugate variables Qn 
and Pn such that: 



and let us write: 



^R,n 



Pn), 
Pn)- 



(C13) 

(C14) 
(C15) 



(C16) 
(C17) 



We have: [$fl„„, $l,„] = ztt - 7r([g„, P„] - [P„, Q„])/2 = 
0. Then the fermion operators are rewritten as: 



'^Ra 



1 



\/27ra 
1 



e --.....-^.e-v 2.-." -.w^^^ (C18) 



(C19) 



(C20) 



Therefore, we can define a set of Majorana fermion op- 
erators: 



-P„) 



fin 



(C21) 
(C22) 



which satisfy the commutation relations (|(]6(l . As a 
result, the representation IjClSp and the representation 
(IC7|) are equivalent to the representation (|C1I) . 

It is well known that one can also define bosonization 
using non-chiral fields^^ (/)„ given by: 



(C23) 

(C24) 



In the case of the bosonization procedure ICip . the 
non-chiral fields have the commutation relation: 



[(j)n{x),(j)yn{x)] = 0, 
[e{x)r,,Om{x')] = 0, 



(C25) 

(C26) 



[9{x)n-,(t>ni{x')] = -i-(5„,mSign(x - a;'), (C27) 

In the case of the bosonization procedure ljC7p . the 
non-chiral fields can be defined similarly but they 
have the commutation relation: [^(a:)^, 0„i(a;')] = 
—i'K5n,mQ{x — x'), whcrc Q is the Heaviside function. 
The relation between the two sets of chiral fields is: 



(t>n{x) = (t>n{x) 



9n{x) = f?„(x) + \/-P„ 



(C28) 
(C29) 



An interpretation of Qn is that it is proportional to 
an operator counting the number of fermions (both right 
and left moving) of type n. P„ is the proportional to the 
phase conjugate to this fermion number, and thus must 
be compactified. 

Note that if we perform a rotation on the fields On 
and (pn, this rotation preserves the commutation rela- 
tions (jC25p . This is a crucial property. It also preserves 
the commutation relations in the case of field theoretic 
bosonization, and obviously it preserves the relation of 
commutation of the zero modes P and Q. Therefore, one 
always goes from the condensed matter convention to the 
field theoretical convention by using the same shift of the 
non-chiral fields (|C28|) by zero modes. 

When converting products of two fermion opera- 
tors from the condensed matter convention to the 
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field-theoretical convention, the following rules apply: 

'nR/L,nVR/L,n' ~^ VnVn' for different species n n' and 
for a given specie n : 



i2(t>„ 



(C30) 



We close this section with two remarks. First, the field 
theoretical representation can be understood naturally 
on a semi-infinite system extending to -|-cxi, with an open 
boundary condition at the origin^^^ : the origin is then 
sent to —oo to recover an infinite system. The non-zero 
value of the commutator of the right and left moving field 
possess then the physical interpretation of left movers be- 
ing reflected into right movers at — oo. As a second re- 



mark let us note that when applying the condensed mat- 
ter bosonization to the double Ising model, one finds that 
the operator cos does not have the same commutation 
with the the fermion field as the products of Ising disor- 
der fields^^, whereas this relation is satisfied in the field 
theoretical bosonization. The reason for this difference 
can be inferred from our first remark : the construction 
of the product of disorder fields as a string is made for 
a system with an open boundary condition at the ori- 
gin. This boundary condition is included from the start 
in field theoretical bosonization but not in the condensed 
matter bosonization. Thus, when using condensed mat- 
ter bosonization, Klein factors must appear in the ex- 
pression of the products of the Ising disorder fields as a 
function of cos 6. 
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